Local-field approach to the interaction of an ultracold dense Bose 

gas with a light field 



Konstantin V. Krutitsky *, Frank Burgbacher and Jiirgen Audretsch * 
Fakultat fur Physik, Universitdt Konstanz, Fach M 674, D-78457 Konstanz, Germany 

(February 1, 2008) 



Abstract 

The propagation of the electromagnetic field of a laser through a dense Bose 
gas is examined and nonlinear operator equations for the motion of the center 
of mass of the atoms are derived. The goal is to present a self-consistent 
set of coupled Maxwell-Bloch equations for atomic and electromagnetic fields 
generalized to include the atomic center-of-mass motion. Two effects are 
considered: The ultracold gas forms a medium for the Maxwell field which 
modifies its propagation properties. Combined herewith is the influence of the 
dipole-dipole interaction between atoms which leads to a density dependent 
shift of the atomic transition frequency. It is expressed in a position dependent 
detuning and is the reason for the nonlinearity. This results in a direct and 
physically transparent way from the quantum field theoretical version of the 
local-field approach to electrodynamics in quantum media. The equations for 
the matter fields are general. Freviously published nonlinear equations are 
obtained as limiting cases. As an atom optical application the scattering of 
a dense beam of a Bose gas is studied in the Raman-Nath regime. The main 
conclusion is that for increasing density of the gas the dipole-dipole interaction 
suppresses or enhances the scattering depending on the sign of the detuning. 
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I. INTRODUCTION 



During the last decade we have witnessed a rapid progress in the field of atom optics [|TJ . 
Examples for experimental success are Bragg scattering of an atomic beam from a standing 
wave p|, reflection of atoms from an evanescent wave 0, and interference of atomic beams 
0]. Theoretically predicted are among other things the transmission and reflection of atoms 
from standing light fields ||, tunneling |J, and also it seems to be possible to focus atomic 
waves with the help of inhomogeneous fields 0. All these results refer to a physical regime 
in which the system can be described effectively as a single particle system. The respective 
quantum equations of motion contain no interaction between particles, which makes them 
linear. 



With the experimental realization of Bose condensates R- U | , atom optics has reached 
a new stage of development. The influence of many body effects on the atomic motion in 
light fields became important, demanding a quantum field theoretical treatment of matter. 
In addition, in the presence of light the polarizable atoms are subject to the resonant dipole- 
dipole interaction which causes the equation of motion to become nonlinear. Perhaps the 
most well known effect of nonlinear wave propagation is the formation of solitons, which have 



also been studied for condensates in the pioneering theoretical work published in [fL2-14| . For 



the shaping of these solitons see [|T3| . In the framework of this new domain of nonlinear atom 
optics too, atom optical devices were studied JlfJ . For beam splitters for condensates based 
on the Bragg scattering of a coherent beam see [|TT|-|TP|]. The underlying theoretical scheme 
has been criticized in p0[ and PH| , mainly with regard to the treatment of the electric field 
inside the atomic sample of polarizable atoms. Treating the dipole-dipole interaction as a 
contact interaction is regarded as a too severe approximation. In addition, an approximation 
restricting to low densities of the condensate is made. 

To overcome these problems, more rigorous approaches to describe the joint dynamics 
of atoms and light in the quantum regime have been developed. First without the inclusion 



of atomic motion [32]. Later, with applications to atom optics in mind p0 |, generalized 
Maxwell-Bloch equations were derived which in principle include retardation and therefore 
correctly describe dispersion, as well as the local-field effect, and accordingly give the correct 
Maxwell equations. In |S| the formalism was used to calculate the refractive index of a dilute 
Bose gas and in [F2TJ atoms moving in an evanescent wave and the gravitational field were 
considered. We want to contribute to this field of research. 

The second subject to which this paper is related are the local-field effects. According to 
the Lorentz-Lorenz relation the local electric field in a medium differs from the macroscopic 
field. This leads to the local-field effect in classical, quantum and nonlinear optics which 
manifests itself mainly if the density of the molecules, which constitute the medium, is high. 
In classical optics the Lorentz-Lorenz relation leads to the more general Clausius-Mossotti 
formula for the refractive index of the medium in comparison with that based on the standart 
response theory. In quantum optics the local-field effect leads to an essential modification 
of the lifetime of the atom embedded in a dense dielectric medium [E1,E9], to statistical 
effects in superfluorescence and amplified spontaneous emission P6| . Many new interest- 



ing phenomena caused by the local-field effect can be found in nonlinear optics, where the 
local-field correction gives rise to the nonlinearity in Bloch equations and as a consequence 
to the optical bistability that is intrinsic to the material and doesn't require an external 
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feedback (see, for example, |27H30| and references therein), to intrinsic self-phase modula- 



tion in self-induced transparency |31|J3^] , to linear and nonlinear spectral shifts |p3| , M , to 
great enhancement of the index of refraction with no absorption in the systems exhibiting 



coherence- induced lasing without inversion [p5|, to ultrafast intrinsic optical switching in a 



dense medium of two-level atoms 0] and to a nonlinear transparency of a thin resonant 
layer remaining near to its ground state ||37|| . 

In the present paper we shall show that the local-field effect plays also an important role 
in the atom optics of ultracold Bose gas. It causes the nonlinearity in the Bloch equations for 
the atomic quantum field. This leads in its turn to the nonlinear phenomena in atom optics. 
Our purpose is twofold: With regard to the development of a general scheme we study in 
a first part an electromagnetic field propagating in an ultracold gas (Bose condensate) of 
neutral polarizable atoms. Our intention is to give a straightforward and short derivation of 
the coupled dynamical equations for the Maxwell field and the matter field operators taking 
into account the atomic center-of-mass motion. The resulting equations form the basis of the 
respective nonlinear atom optics in the domain, where the resonant dipole-dipole interaction 
is the dominant effect. Starting from the Hamiltonian in second quantization we present a 
physically intuitive scheme which reflects the fact that we are dealing with electrodynamics 
in media. We are thereby directly led to the related concept of the local-field approach. 
The Lorentz-Lorenz formula and the Clausius-Mossotti relation follow immediately, so that 
the optical properties can be discussed. The origin of the nonlinearity of the matter field 
equation becomes transparent. Important is also, that no approximation with regard to the 
density will be made. Previously considered equations can be obtained as limiting cases for 
low density. It is not necessary to refer to a particular quantum state of the gas. 

The second purpose of the paper is to apply the general results obtained for the dynamics 
of the field operators to a particular device. In principle we have to solve a coupled system of 
equations. Under certain assumptions it is possible to solve the Maxwell equations without 
restricting the density of the gas. We study the scattering of the atoms from a standing light 
wave in the strong coupling regime. The dependence on the density of the gas is discussed. 



II. THE MODEL 

We discuss the combined system of a Bose gas of identical neutral polarizable atoms 
interacting with an incident monochromatic field Ej n (r, t) with frequency ujl. The probe 
field E in is a classical external field. No operator is related with it. It is fixed and therefore 
not subject to a dynamics. For example, it may be a standing or running laser field. The 
quantized center-of-mass motion of the atoms is taken into account. The atoms are repre- 
sented as point-like neutral two-level atoms with states |1) and |2) with "bare" transition 
frequency u a . The atoms interact via the resonant dipole-dipole interaction of the induced 
point dipoles. This interaction is mediated by the exchange of photons with annihilation 
and creation operators CkA and cj^ referring in lowest order to plane waves in vacuum with 
wave vector k and polarization A. 

The total atom-field Hamiltonian may then be written in the electric dipole approxima- 



tion in the length gauge or dE coupling 38 in the form 



H = H A + H F + Hai + H AF 
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h a =y: 



h 2 v 2 



+ ?iUJ a (Tjcri } , 



2m 

dJ2 (cri + crj) E in (r h t) 



(1) 



kA 



H 



Al 



Haf = -ft 9k\cl\ ex P Hkr) (a { + <rj) + H.c. , 



i kA 



i^A is the Hamiltonian of the freely moving atoms, Hp is the Hamiltonian of the free photons 
and Haj describes the interaction of the atoms with the incident field E in . The transition 
operator is = |2)(1|. The exchange of photons between the ground and excited states of 
the atoms is incorporated by the atom- field interaction Haf- A rotating wave approximation 
has not been made, d = d* is the matrix element of the atomic dipole moment, g^x 
represents therefore the coupling strength of an atom to the photons. Transitions between 
the internal states of the atoms go back to the coupling of the atom to the incident field on 
one hand, and to the emission and absorption of photons on the other. 

We assume that all contact interactions can be neglected. In order to find the conditions 
under which this approximation is valid, it is necessary to estimate the ratio Ud/U g , where Ud 
and U g are the mean energies of the dipole-dipole interaction and the ground-state collisions, 
respectively. Treating the ground-state collisions in terms of s-wave scattering and restricting 
to dipole optical transitions we get the following inequality |21] 



Ud 3 s 
U n 



8 a± n 



(2) 



where s is the saturation parameter of the atomic transition, a s is the scattering length and 
k a — w„/ c. The scattering length a s is a quantity of the order of 1 nm. In the optical domain 
k a is of the order of 0.01 nm" 1 , and we get the following estimate 



Ud 

u n 



> 37.5s 



(3) 



Thus, the ground-state collisions are negligible if the saturation parameter s is of the order 
of 0.01 or higher. This shows, that we don't need very high values of the parameter s in 
order to neglect contact interactions. 

To take into account quantum statistical effects for the two-level atoms, we turn to 
quantized matter fields represented by two-component vectors 



V(r,t)=Vi(r,t)|l> + V2(r,t)|2) 
The respective bosonic equal time commutators are 



(4) 



o 



SijS (r 



(r-r'), z,j = l,2. 



(5) 
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In second quantization the Hamiltonian operator (|l|) takes then the form: 

2 , / 



2m 



4'(r) + /" rfr^(r)^ a ^ 2 (r) + # F 



(6) 



- J drdE m ^(r)^ 2 (r) - K J dr^^c^exp (-ikr) {^(r)^(r) + $(r)A(r)} - if. 

kA 



C. 



Making use of the Hamiltonian (^) and the commutator relations we obtain the following 
Heisenberg equations of motion for the atomic and field operators: 



ih 



dt 



2m 

n Y,9k\ci\(t)exp(-ikr)^ 2 (r,t) - hip 2 (r, t) ^ fl'kA exp (ikr) c kA (t) , (7a) 



kA 
2V72 



kA 



ih 



difofat) h 2 W 



dt 



4 2 (r, t) + hujj 2 (r, t) - dE+ (r, t)^(r, t) 



2m 

^Hfi , kAC k AW ex P(- ikr )^i( r ' t ) -?#i(r,t) 9k\ exp (ikr) c kA (t) , (7b) 



kA 



kA 



i/i 



dc kA (t) 

<9t 



huj k c kX (t) - hgl x / (irexp (-ikr) ^(r, t)^i(r, t) + tpl(r,t)ip 2 (v,t) , (7c) 



where are the positive and negative frequency parts of the incident classical electric field. 



III. THE LOCAL-FIELD APPROACH 

From the physical point of view we are in fact investigating the propagation of an inci- 
dent electric field through a particular medium. It is therefore reasonable to try to recover 
the fundamental quantities and relations of the well-known theory of electromagnetic fields 
in media [j39,i0|. There one regards the medium as being composed of interacting induced 



point dipoles in otherwise empty space. The incident field E in (r, t) causes the dipoles to 
emit radiation. The electric field to which a single dipole at the position (r, t) effectively 
reacts is then the resulting microscopic local field E/ oc (r, t). On the other hand, the Maxwell 
equations in media refer to a macroscopic or mean electric field E mac (r, t) which is obtained 
by averaging over a region which contains a large number of atoms. E mac (r, t) and E; oc (r, t) 
differ by the local-field correction or Lorentz-Lorenz correction. This correction is propor- 
tional to the macroscopic electric polarization P(r, t) which also appears in the Maxwell 
equations. We will now show that the dynamical operator equations which we obtain from 
(0) reflect in fact the physical structure of this local-field approach. 
The formal solution of for the photon operators is 

PkA(t) = |ckA(0) 



rt 



10 



+igtx I dt' I dr'exp (iu k t' - ikr) ^|(r', t')^^', t') + $(r', t')ip 2 (r', t') \ exp (-iu k t) , (8 
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where the first term CkA(O) refers to the free-space photon field and the second one goes back 
to the interaction with the atoms. 

To study the back reaction of the photons on matter we insert fl8|) in ( |7a|) and ([T5|) and 
rearrange the matter operators with the help of the commutator relations @. For details 
of the calculations see Appendix A. Introducing the polarization operator P^r, t) of the 
ultracold ensemble 

P ( r, t) = d$(r, t)$ 2 (r, t) + H.c. = P+(r, *) + P~(r, *) , (9) 

we obtain as a central result the dynamical equations for the operators of the two matter 
fields 

~ h Yl 9tx d lx(°) ex P (~^ kr + iuJ kt) ^ 2 (r, t) 

kA 

- h^ 2 {r,t) ^fl'kACkA(O) exp (zkr - iu k t) , (10a) 

kA 

ih ?^?A = -^Mr, t) + h{u a + 5- vy/2) Mr, t) - dE+ c (r, t)^(r, t) 
9k\el\( Q ) ex P (-^ kr + iuJ kt) $1 (r, t) 

kA 

- hi>i(r,t)J29k\Ck\{0) exp (ikr - zo; fe t) , (10b) 

kA 

where we have introduced 

E£c(r, = E^(r, t) + / dr'V x V x P± (r ^~ R/c) , (11) 

and Vx refers to the point r. 

Eq. ([LTD shows that E toc (r, t) is the superposition of the incident field E in (r, t) and the 



electric field which is obtained by adding up at the space-time point (r, t) the electric dipole 
radiation of all other atoms. E^ c is therefore the operator of the local field. It is correct for 
all distances including the near-field zone. As one should expect it is this local field which 



drives the inner atomic transition in (10). Because of this fact the matter field equations 
(PH|) together with (JTTI) and (|9|) are nonlinear. This follows directly without approximation 
from the Heisenberg equations (^). The physical reason is that the electromagnetic field 
which propagates through the medium causes the atoms to emit and absorb photons of the 
dipole radiation. 

Note that in the eq. ([LTD a small volume around the observation point r is excluded from 
the integration. The eq. (|TT|) can be written down in another form, where the integration 



is carried out over the whole space, besides there is an additional term in the r.h.s. which 



containes a 5-function |23| , |T0"|| . These two forms of integro-differential equation are totally 
equivalent and give the same results. 

The equations ( |T0"D are nothing but the atom-optical analogue of the usual optical Bloch 



equations |41,22|. They describe the dynamical evolution of the second quantized matter in 



the field of electromagnetic radiation. In the Appendix B we shall show that under certain 
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conditions the optical Bloch equations which describe the dynamical evolution of the first 



quantized matter can be derived from the equations (|10|). 

Treating the quantum g medium, it is important to establish a connection to 

the observables of the quantum electrodynamics in media. Having an application to an 
ultracold gas in mind, we have to relate the local field E; oc (r,t) to the macroscopic (or 
mean) field E mac (r,t) which is obtained by averaging in space over a region which contains 
a great number of atoms. For a very low density ultracold Bose gas the region on which the 
macroscopic field is averaged may contain actually fewer than one atom, yet the averaging 
technique should still be valid due to atomic derealization. 



We generalize the idea of Lorentz as it is described in p^ , |39| , P0|] and assume the local 
field to be a sum of several contributions: 

E ioc (r, t) = E near (r, t) + E mac (r, t) - E P (r, t) . (12) 



■'near 



Let's consider a macroscopically small, but microscopically large volume V around r. E r 
is a contribution from the dipoles within V. E mac is the macroscopic field obtained from the 
field of all dipoles in the medium. Because the dipoles in V are already taken into account in 
E n ear, we have to subtract from E mac the contribution Ep of these dipoles as it is obtained 
in the averaged continuum approximation. Ep is therefore a macroscopic quantity. It can 
be related to the polarization, which is an averaged quantity too, according to 

47T - 

Ep(r,t) = -— P(r,t). (13) 



For a novel derivation see [p2]| . In the following we restrict to the purely classical local-field 



correction which corresponds to a vanishing E near (comp. p4p9|,fi0| ). Thus, we see that, 
because of the influence of the near dipole-dipole interactions, the local field is obtained 
from the macroscopic field in adding the local-field correction (comp. e.g. [ p^ , j39| ,[4C|1). 



4-7T 

Ef oc (r,t) = E± ac (r,t) + y P ± (r,t). (14) 

This equation is also often called the Lorentz-Lorenz relation. The propagation of the op- 
erator of the macroscopic field is thereby given by Maxwell wave equation for a charge-free 
and current-free polarization medium 

VxVxE^ (rt)- i gg^M) 47r ffP*(r,t) (15) 
Vx VxE mac (r,t)- ^ m2 ^ m2 . (15) 

When taking the expectation value, the operator P(r, t) leads to the averaged dipole moment 
per unit volume. 



To sum up: Taking together eqs. (fL5|) and the truncated eq. (0) with (|14|) , we have 
obtained a self-consistent system of coupled differential equations for the collective dynamics 
of the fields E mac (r,t) and ^1,2(1*, t). Note that when dealing with the Maxwell propagation 
equation (|15|), the existence of the incident field E in has to be taken into account. It 
initiates the polarization of the medium. With ([14]) and (0) the matter field equations (10) 



are nonlinear. This is the direct consequence of the fact that matter couples necessarily to 
the local field as the driving field for internal atomic transitions and not to the macroscopic 
electric field. 
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IV. NONLINEAR MATTER EQUATION AND OPTICAL PROPERTIES 



A. General case 



In the next step we decouple the dynamical equations for the two matter field operators. 
We substitute (0) in ( |lUa| ) and ( |10bp and pass to the reference frame rotating with the 
frequency ujl of the incident field which is assumed to be monochromatic 



to obtain 



ih 



dt 



th- 



at 



E„ ac (r, t) = £+ ac (r) exp (-iu L t) , 
ip 2 (r, t) = 2 (r, t) exp (-iu L t) , 



2m 2 3 

-^0 2 - 5n + (r)^ - f dHXU* 



-h (A + i 7 /2) 02 + G 2 



(16) 



(17a) 



(17b) 



with the detuning A = u;^ — uj a — 5. The position dependent Rabi frequency Q (r) = 
2d£^ ac (r) /h is related to the macroscopic electric field. The noise terms in the eqs.(|T7j) are 



G 1 (r,t) =-^[fi(r,t)0 2 + 2 f 2 (r,t) 
G 2 (r, t) = -h fj (r, t) ^ + i (r, t) 



(18) 



with operators 



fi (r, t) = fi'kACkA(O) exp [-i (u k -u L )t + ikr] 

kA 

f 2 (r, t) = 9k\c kX (0) exp [-% (u k + uj L )t + ikr] 



(19) 



kA 



which give the effect of vacuum fluctuations on the atomic quantum field. 

We assume sufficiently large detuning such that the spontaneous emission gives a small 
contribution. We may therefore apply the adiabatic approximation to ( 17bD (comp. |Ti 



2 (r,t) 



n+(r)Mr,t) 



G 2 (r,t) 



2 (A, (r, t) + i 7 /2) h (A, (r, t) + * 7 /2 



(20) 



and eliminate the field ip 2 from our scheme. The position dependent local detuning introduced 
in© 



47T 

A,(r,*) = A + — dVI(r,t)Vi(r,t) 



(21) 



is an operator depending on the density operator of the ground state. We interpret this as 
a shift of the internal energy levels which may increase or decrease with increasing density 
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depending on the sign of the detuning A. Corresponding frequency shifts can be found in 
nonlinear optics where they lead to the effects mentioned in the introduction. Substituting 
(pl|) in ( |17a| ) we obtain as the intended result the nonlinear dynamical equation for the 
matter field operator tpi (r, t) 



ih 



dt 



h 2 v 2 h 
+ 



A - 27/2 



2m 4A i 2 (r,t)+f/4 



v R = h- 



A-Z7/2 
A?(r,t)+r74 



I(r 2 + rt)^ + i 



1 2 



r 1 )n-+ ri + r 2 



(22) 



It couples to the macroscopic electric field via the Rabi frequency f2 ± (r). Here Vr is a 
random potential caused by vacuum fluctuations. Its vacuum average (Vr) vanishes. 

On the assumptions stated in Sect. I, the equation fl22|) is general. No approximation 
apart from the adiabatic approximation has been made. The dipole-dipole interaction me- 
diated through the exchange of dynamically produces photons is completely contained. The 
influence of the atoms is taken into account in the Maxwell equations fll5|). The local-field 
correction is included. It is not assumed that the density is low. The nonlinearity goes 
back to the density dependent local detuning Ai(r,t) of eq. fl2~Tl). For increasing density 
and positive detuning A the local detuning grows and correspondingly the nonlinear term 



in (|22|) representing the coupling to the macroscopic electric field becomes smaller. On the 
other hand, for negative detuning the absolute value of the local detuning decreases and the 
nonlinearity becomes greater. When (p2|) and the Maxwell wave equation fll5|) are solved, 
all possible effects on the center-of-mass motion of the atoms can be worked out. 

Because we are mainly interested in atom optical problems and want to study the coher- 
ent evolution of the center-of-mass motion of the gas, we shall neglect spontaneous emission. 
This is valid for situations where the absolute value of the local detuning is much bigger 
than the spontaneous emission rate |A/| 3> 7. In order to do this we drop in the following 
the random potential Vr and the spontaneous emission rate 7 from our equations. 



B. Special cases 



The equations which are usually used for the description of an ultracold ensemble put in 
an external laser field may be derived from (|2~2~D as limiting cases. Neglecting dipole-dipole 
interactions we have Aj = A and E^ c = E^ ac = E^. The eq.(p2[) reduces to the well-known 



equation of the single particle theory [42 



ih 



dt 



h 2 w 2 



h 

2m 4 



&(r,*) 



(23) 



Expanding ( p2|) to the lowest order in the density tp\ipi we get an equation of the Gross- 
Pitaevskii type which is here obtained in the local-field approach: 



ih 



dt 



h 2 v 2 

2m 



A 



2nd 2 



$(r,t)&(r,t) 



A(r,t) . 



(24) 
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It may be applied fow low densities. A similar equation has been obtained in [17 



Another approximative approach is the following: Introducing the standard (linear) 
atomic polarizability a = —d 2 /hA, we can rewrite the nonlinear term in eq.(p2|) as 



h 
T 





6+(r) 


2 


A 


1 - ^ot$$ 





(25) 



Neglecting the term in the denominator proportional to the density squared, we get 

2 

Mr,t), (26) 



ih 



Of 



h 2 v 2 h 
+ - 



2m 



4A 1 



which agrees with the eq.(3.12) of the paper [21 for vanishing contact interaction and 
vanishing gravitational potential apart from a factor 2 (we have 87r/3 instead of 4tt/3). 



C. Optical properties 

We return to the general case of Sect. IV to discuss the optical properties of our ultracold 
Bose gas. We insert the adiabatic solution (p0|) into the polarization (|9|) to obtain 

P + M) = xM)E+ QC (r,t) (27) 

and P~ correspondingly with the operator of dielectric susceptibility of our medium 



q^|(r,t)^i(r,t) 
l-fa$(r,t)^(r,t) 



where a is the atomic polarizability stated above. ( p7|) with (|2q ) corresponds the well-known 
Lorentz-Lorenz formula. Inserting (^71) in (JTS[) we find together with ( |2"2"| ) the new system of 
coupled equations for the quantum fields E^ ac (r, t) and ipi(r, t) as it is valid in the adiabatic 
approximation. 

Assuming for this system that the spatial variations of the atomic density are not very 
large so that div'Ei^ nac ~ 0, we have 

V 2 ^ ac + k 2 L h 2 £± ac = , (29) 

whereby the refractive index satisfies the Clausius-Mossotti relation: 



1 + tfa$lij) 
1 - 4ra$^i 



n 2 = . L r V (30) 



Based on our system of coupled equations and the related optical interpretation, we turn 
now to an application and discuss an atom optical device for ultracold atoms, which can be 
used as a simple model for a beam splitter. 
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V. DIFFRACTION OF AN ULTRACOLD ATOMIC BEAM FROM A STRONG 

STANDING LIGHT WAVE 



We have to solve both the equation of the quantum field (|22|) and Maxwell's equation 
(f29l) with (|30|) . In general this has to be done in a self-consistent way and it will be difficult 
to do so analytically. Nevertheless, some simplifying approximations can be made when 
applying the whole formalism to the diffraction of an ultracold atomic beam by a standing 
electromagnetic wave. For this we generalize and modify the approach of JlT] . 



We consider the following setup: An incident atomic beam moves in ^-direction, per- 
pendicular to a standing laser beam consisting of two counter propagating waves along the 
y-axis with wave vectors +nk.i and — nh^, respectively, and with a Gaussian envelope. The 
motion in ^-direction remains classical during the whole evolution and only the diffraction 
in y-direction is treated quantum mechanically. In order to get a distinct diffraction pattern, 
we assume that the width of the atomic wave packet w y is sufficiently large. Additionally 
the variations of the density should be small on a scale of the effective wavelength of the 
laser (w y ^> 2tt /nk L ). In this case the atoms can be described as a homogeneous medium. 
The effect of the atoms on the laser beam is purely dispersive and only the wavelength will 
be shifted. The solution to ( [29D with (|30D is then given by 

Vt + = \{l \ 2 exp (^—z 2 /w 2 ^) cos 2 nk L y . (31) 

To solve the equation of motion for the atomic quantum field, we assume that we are in 
the Raman-Nath regime and that in equation fl22|) the kinetic term can be neglected during 
the interaction of the atoms with the electromagnetic field. This approximation is valid for 
heavy atoms or if the interaction is so strong that atoms can take up momentum without 
changing considerably the velocity. In this case the density of the atoms remains unaltered 
and their phase changes. With the definition 

V,4^ (32) 
3ft A ' K ' 

of the characteristic volume, we then get for the equation of motion of the atomic beam 



Mi 1 



i- 



dt 4A (1 + V oPl ) 



k , (33) 



where p\ = ipl'tpi is the density of atoms in the ground state. This also guarantees that 
the conditions leading to fl3i~l) are always satisfied, because the density distribution remains 
constant in time and only the phase of the atomic quantum field changes. Making use of 
the assumption that the motion in z-direction is purely classical, we change the variable 
t = z/vg in (|33| ) with v g being the group velocity of the atomic beam. Then the solution of 
eq. (|33"D for z ^> Wl (in the far zone) may be written in the following form 

( f°° -i l^ + | 2 \ 
^!(y,oo) =Vi(y,-oo)exp / — — — — ^dz . (34) 

\J-co AAv g (1 + V pi) J 
Representing pi as a Gaussian wave packet with width w y 
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we substitute ([Jl]) and fl35|) into (|34] ) and get after integration 



$i(y, oo) = ^i(y, -oo) exp 



pi = p exp (-y 2 /wf) , (35) 



(36) 



— i4g cos n/cL?/ 



1 + V po exp (-y 2 /^ 



where g = -§£^V^- 

The relative broadness of the wave packet compared to the laser wavelength can be 
expressed as w y ^> 2ir /nkz,. Herewith we can represent the solution (^) in the form of a 
Fourier series expansion: 

oo 

My, oo) = ^ (y, -oo) ]T e i2 ^ fc ^(-i)V ? (r) , (37) 

q= — oo 

where r = 2g / (1 + Vopo) 2 and Jq is the g-th order Bessel function. Accordingly the prob- 
ability P q to find the wave packet in a momentum state qnk^ is given by 

p q = J q( T )i Q = -oo,...,+oo , (38) 

with Pq corresponding to the intensity of the outgoing beam which propagates in the same 
direction as the incident atomic beam. The angle of diffraction a q for a particular momentum 
state is thereby given by 

qnhk,L /0 „x 
tana„ = . (39) 

ITLVg 

Therefore the diffraction pattern as it follows from (|3~8"D depends on the initial density po of 
the beam which is contained in r. If po = 0, the formula ( |3"8"| ) coincides with the formula 
(14) of paper 0. 

The result (|||) provides the possibility to obtain semi-quantitative estimates of the 
influence of the initial density po of the ultracold atomic beam on the intensities of the 
various diffraction orders. Let's consider first the case of positive detuning A. It is obvious 
that in this case if po increases, the dimensionless effective interaction time r will decrease. 
In addition, taking into account the form of the Bessel functions, the relative intensities of 
the diffraction orders with q ^ decrease and Po increases. In the case of negative detuning 
A the absolute value of r decreases with the increase of the density p . This implies that the 
effect of increasing density in the case of negative detuning A is to increase the coupling to 
the laser, enhancing the diffraction pattern. Physically these features originate from a shift 
of the internal energies of the atoms going back to the interaction with the other nearby 
atoms, as it is clearly expressed by equation (|2T|). The detuning of the laser beam changes 
locally with the density. With increasing density the coupling to the laser in equation 
(f22|) diminishes or increases depending on the detuning A. Thus, our estimates show that 
with the increase of initial density, the atomic beam wave properties being manifested in 
the diffraction can be suppressed or enhanced due to dynamical dipole-dipole interactions. 
These effects are significant when Vopo ~ 1, i-e. when po ~ (A/^k^/ir. In the optical 
domain k\ is of the order of 10 15 cm~ 3 . Taking into account that due to the fact that in 
our analysis | A| ^> 7 the estimated value of the density required for large nonlinear effects 
at present is clearly larger than what can be obtained for alkali atoms in a magnetic trap. 
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VI. CONCLUSION 



We have treated the interaction of a dense Bose gas with light in the local-field approach. 
As an application we have shown that it is possible to obtain analytical results for the the 
problem of scattering of an ultracold gas from a strong standing laser wave in the Raman- 
Nath regime. These results are not limited to low densities and include the dipole-dipole 
interaction consistently. Two effects had to be considered: First, the dynamical dipole- 
dipole interaction of nearby atoms shifts the atomic transition frequency. This gives a 
locally defined detuning between the atomic transition and laser frequency, which depends 
on the local density. Secondly, the atomic gas acts as a medium for the electromagnetic 
field and has to be taken into account in the Maxwell equations. While retardation of the 
electromagnetic field can always be neglected for short distances and gives only a negligible 
contribution to the level shift, retardation plays an important role at higher densities, and 
is necessary to describe the dispersion of electromagnetic waves. 

Both effects give a contribution to the center-of-mass motion of the atoms in the gas 
and have to be considered simultaneously if the formalism is applied to atomic scattering. 
If variations of the density of the gas are only on length scales which are long compared 
to the laser frequency, the refractive index of the gas can be defined. The corresponding 
expression coincides with the well-known Clausius-Mossotti relation, and allows oneself to 
find a solution to the Maxwell equation alone, provided the conditions for the validity of 
the solution are preserved during the scattering process. This is guaranteed in the Raman- 
Nath regime, where the motion of the atoms is neglected and therefore the density remains 
constant. The result is again a standing wave, but the wavelength has changed according to 
the refractive index. With this, the recoil on the atoms is different compared to the case of 
no medium considered. Additionally, because of the density dependence of the atom-laser 
detuning, the time scale of scattering may be prolonged or decreased for high density gases, 
giving a suppression or enhancement of the scattering. 



ACKNOWLEDGMENTS 

This work has been supported by the Deutsche Forschungsgemeinschaft and the Op- 
tikzentrum Konstanz. One of us (K.V.K.) would like to thank the members of the AG 
Audretsch at the University of Konstanz for many interesting discussions and kind hospi- 
tality. 



APPENDIX A 



To derive the eqs. (|i0|) and (|TT|) which are a central result of the local-field approach, we 
insert the solution (||) for photon operators into the eqs.(^) for the matter field operators. 



From (7a) we obtain for the field operator of the lower state 



dipx(r,t) 



2V72 



rv 



■Vi(r,t)-dEi;(r,t)^(r,t) 



ih 

at 2m 

9k\cL\(fy exp (-ikr + iu k t) Tp 2 (r, t) - /# 2 (r, t) 9k\ exp (ikr - iu k t) c kA (0) 



kA 



kA 
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+X(r,t) + Y(r,t)+Z(r,t) 



(40) 



where 



X(r, t) = -ih I dr' f dt'^\(r', t')Mr' , f) £ \g^\ 2 {exp [ik (r - r') - iu k (t - t')} - c.c.} 



xip 2 (r,t) 



(41) 



Y{r,t) = -ih / dr' 



dt' £ |s kA | 2 exp [ik (r - r') - iu k (t - 0] Mr, t)4(r' } t')Mr', t') 



kA 



(42) 



Z(r, t) = ihf dr' f dt' I^aI 2 exp [-ik (r - r') + iu k (t - t')} $(r', t')Mr', t')$ 2 {r, t) . 

(43) 

We turn to the evaluation of X(r,t) first. With 



(44) 



and 



(de kA ) 2 = £ ^eg^e^ 



(45) 



m,n=l 



we find 



X(v.t) = - l -jdr'^dt , $ l {r',t')i) 2 {r'it') £ d m d r , 



m,n=l 



X 



27fflU) k 



E exp [ik (r - r') - iu; fc (t - t')} - c.c. \ ^ 2 (r, f) . (46) 



kA 



Introducing the commutators of the spatial components E m (r, t), m — 1, 2, 3, of the operator 
of the free electric field (comp. section Cm. 2 of f3"8||), we may write: 



X(r,t) = ~ [ dr' fdt'^lir', t')^ 2 (r' ,t') £ d m d n [E m (r, t), E n (r', t')] 4> 2 (r,t) . (47) 



As can be found in the section Cm.3 of @] the commutators may be rewritten in the 
following form: 



E m (r,t),E n (r',t' 



) {3R m R n 



5'{R - ct) - 5'{R + ct) S(R - ct) - 5{R + or) 



R 2 



R 3 



R m R n _ x \S"(R - ct) - 5"(R + CT) 



R 2 



R 



(48) 
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where R = r — r' and t = t — t' . We introduce the polarization operator 

P(r, t) = d$(r, t)$ 2 {r, t) + H.c. = P+(r, t) + P"(r, t) , (49) 
and sum in (^7|) over m and n making use of (f48l) . After integration over time we get: 



X(r,t) = -d / dr' 



3( 



n n 



R 3 






r " +" 




r - +" 




r ~ +i 


n — 


P 




P 


n^j n — 


P 



1? 



R 2 



Mr,*) 



(50) 



with n = R/i? and P + = P + (r', t — R/c). This agrees formally with an integral over the 
retarded fields of electric dipoles situated at (r',t — R/c) (comp. section 2.2.3 of p5[| ) and 
may therefore be rewritten as 



X(r,t) = -d / dr'V x V x -± ^(M) 



(51) 



where Vx refers to the point r. 

To work out Y(r,t) we want to introduce [P~] in the same way as we have introduced 
[P + ] in X(r, t). To do so we have to change the order of the matter fields. With the help of 
a parameter e <C t and with reference to the time r = t — t' we obtain Y(r, t) in the form of 
two integrals: 



Y(r,t) = -ih I dr' £ dr ^ Ig^xf exp [zk (r - r') - iu k r} ip 2 (r, t)fy(r',t - r)^i(r',t - r 



-ih dv' drJ2\g kX \ exp [ik (r - r') - iuj k r] V4(r', t - T)ipi(r', t - r)i(j 2 (r, t) . (52) 

J Je kA 

Making use of the equal time commutators of the matter field ip 2 at r — > we get 
Y (r, f) = -ffl / * f dr £ | JkA | 2 exp [ik (r - r<) - i. k r] t - rtf, (r<, f - r)*(r, f) 



kA 



kA 



-^XH^aI 2 / rfrexp {-iu k r) $i{r,t - r) 



(53) 



We now turn to the sum Y(r,t) + Z(r,t) and repeat the procedure we have described 
above when working out X(r,t). This leads us to a term corresponding to fl51"|) with 

replaced by P~ and a rest. Collecting all terms we then get 



X(r, t) + Y(r, t) + Z(r, t) = -d j dr'V x V x ^^^(r, t) 

-d / dr'V x V x ^— ^V2(r,t) -i^y]bkA| / drexp(-ia; fc r) V»i(r,t-r) . (54) 
7 i2 k A 70 
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We pass to the reference frame rotating with the frequency of the atomic transition u a , 
so that ip2 — 02exp (— iu a i). This shows that the term with [P + ] may be neglected in the 
rotating wave approximation. For t — > oo, e — > oo with e <C t the last integral in (1541) gives 



\9k\\ 



kA 



TlS (uJ k ) — iP 



1 

U k 



Vi(r,t) . 



(55) 



Calculating this term in the limit V — > oo for the quantization volume we see that it 
vanishes in accordance with the fact that there is no decay rate an Lamb shift for the lower 
atomic level. Hence ([53]) and ( f4"U| ) finally give the intended result ( |10a ) after combining the 
polarization term with the incident field E^(r, t) to the local field E ; ~ c (r,t) according to 
(0). 

For the matter field ip2 of the upper state we have to repeat the calculations starting 
with (|7aD and (H) to obtain Eif oc (r,t), correspondingly. After going to the rotating reference 
frame, the last integral on the r.h.s. of (1541) becomes in this case 



kA 



2 f°° 

-ihJ2\gkx\ / drexp [-i (u k - u a ) r] </> 2 (r,t) exp (-iu a t) 
kA Jo 



ih^2 \9k\\ 



kA 



tt5 (u k - u a ) - iV 



u k — u a 



which leads to the Lamb shift S and decay rate 7 in (|10b|). 



(56) 



APPENDIX B. THE DERIVATION OF OPTICAL BLOCH EQUATIONS 

In the present paper we are developing the quantum theory of light-matter interaction 
for second-quantized matter. It is possible to show, that the equations which are used in the 
usual nonlinear optics (where the matter is first-quantized) can be derived from our general 
equations fllPl). 

Neglecting the center-of-mass motion and noise terms, we rewrite the equations ( |TUD in 
the following manner: 

= _d ^ c(r ' t) ^ (r ' t] ' (57a) 

^dt^ = h ^ + 6 - *T/2) A(r, t) - dE+ c (r, t)^(r, t) . (57b) 

In the reference frame rotating with the frequency uj^ of the incident electric field the 
field operators have the form 

E L = £?oc ex P (~ iuJ Lt) , ^2 = 02 exp {-iu L t) , (58) 
and we obtain for (|57D 
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ih 



at 



-d£f oc (r,t)4> 2 (r,t) 



ih- 



Of 



-h (A + 2 7 /2) ^(r, t) - d£+ c (r, t)^(r, t) 



(59a) 
(59b) 



with the detuning A = uj l — uj a + 5. 

In the rotating reference frame the macroscopic polarization operator has the form 

P+ = dpR + exp (-iuJ L t) , (60) 

where pR + = iplfo and p = VWi + 0202 is the total density of the atoms which is a sum of 
the densities of the atoms in the ground and the excited state. 

Making use of the eqs. ( |59|) one can show that the time evolution of the operator R + is 
described by the following equation 



^ = (zA - 7 /2) R+ - Z -d£+ c W 



(61) 



where 

w = (j>U 2 - 4>Ui) IP (62) 
is the inversion operator. Its time evolution is governed by the Heisenberg equation: 



dW 

~~dT 



7. 



1 + W) 



2i 
J 



dR + £, 



loc 



(63) 



In the derivation of the equations (|6lD , (|63| ) we have taken into account that the total 
density of the atoms p is an integral of motion. It is possible to show this starting with 
the second-quantized Hamiltonian (Q) where the first term, which describes the center-of- 
mass motion, is neglected. However, it is not possible to show this when starting with the 



equations (|57]) because in these equations the vacuum fluctuations of the photon field are 
neglected. The neglect of the vacuum fluctuations in the equations fl57D corresponds to the 
neglect of some terms which are not presented explicitly in the Hamiltonian ([]). 

Replacing the spontaneous emission rate 7/2 in the equations (|6TD , ( |63|) by the phe- 
nomenological damping and relaxation terms, ji and jt, respectively, we get the optical 
Bloch equations 
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